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Ramin G. Daghigh^ and Yutaka Hosotani^ 

^Department of Applied Mathematics & Theoretical Physics, The Queen's 
University of Belfast, Belfast, BT7 INN, Northern Ireland, UK 
^ Department of Physics, Osaka University, Toyonaka, Osaka 560-0043, Japan 



We investigate gravitating lumps with a false vacuum core surrounded by the true vacuum 
in a scalar field potential. Such configurations become possible in Einstein gravity in the 
presence of fermions at the core. Gravitational interactions as well as Yukawa interactions 
are essential for such lumps to exist. The mass and size of gravitating lumps sensitively 
' depend on the scale characterizing the scalar field potential and the density of fermions. 

These objects can exist in the universe at various scales. 

^ ■ §1. Introduction 

m ■ 

^ I When a scalar field potential has two non-degenerate minima, the absolute min- 

imum of the potential corresponds to the true vacuum, while the other to the false 
C3 ' vacuum. Systems in which false and true vacua coexist are of great interest. The 

I universe is full of various structures such as black holes, stars and dark matter. Does 

• there exist any structure that has a false vacuum in its core? If the entire universe 

O ■ is in a false vacuum, it decays into the true vacuum through bubble creation by 

"q ' quantum tunneling.^) What happens if a false vacuum core is surrounded by the 

true vacuum?^' ' If the size of the core is smaller than the critical radius, the core 
5—1 . would quickly decay, with its energy dissipating to the spatial infinity. If the size of 

?f^| the core is larger than the critical radius, the core becomes a black hole. In either 

_ ^ ■ case the configuration cannot be static. The fate of false- vacuum bubbles has been 

^ , extensively discussed in the literature. 

^ I It has been shown recently that new configurations, cosmic shells, emerge in 

- - ' a simple real scalar field theory in which spherical shells of the true vacuum are 

immersed in the false vacuum."^) Such static cosmic shells can exist thanks to grav- 
itational interactions. However, a static ball of the false vacuum immersed in the 
true vacuum is not possible. 

In this paper we demonstrate that a static false vacuum core becomes possible if 
there is additional matter coupled to the scalar field such as fermions. The number 
of fermions necessary is not large. Gravitational interactions as well as Yukawa 
interactions play a key role in making such a structure possible. It is a gravitating 
fermionic lump. We stress that gravitating lumps are quite different from Q-balls, 
boson stars and Fermi-balls. In Q-balls the conserved charge of the scalar field plays 
a central role in realizing the stability.^-* There is no such charge of the scalar field in 
our model. In boson stars, gravitational interactions as well as the conserved charge 
play a key role.^^'^*^^ The model for Fermi balls is very similar to ours^^^'^^^ in which 
the false vacuum is surrounded by the true vacuum. In Fermi balls, fermions are 
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localized in the transition region, or the domain wall, i.e. they reside in the surface 
region of the lumps. In our model fermions reside in the bulk region inside the lump 
and the Yukawa interaction becomes essential. 

Gravitational interactions can produce lump solutions in non-Abelian gauge the- 
ory as well. Even in the pure Einstein- Yang-Mills theory stable monopole solutions 
appear in the asymptotically anti-de Sitter space. ^^^"^^^ We show in the present 
paper that such gravitational lumps appear even in a simple scalar field theory with 
fermions. 

The spacetime geometry of gravitating lumps is either anti-de Sitter-Schwarzschild 
or de Sitter-Schwarzschild. Dymnikova has discussed the global structure of de Sitter- 
Schwarzschild spacetime, assuming an appropriate external matter distribution.^^) 
The existence of gravitating lumps in the present paper shows that anti-de Sitter- 
Schwarzschild spacetime is indeed realized in a very simple system. Its existence can 
be inferred from the energetics consideration in flat spacetime as well. The exten- 
sion to de Sitter-Schwarzschild spacetime is reserved for future investigation. One 
may explore such objects in the universe at various scales and epochs. We shall see 
that the size of gravitating lumps sensitively depends on the scale characterizing the 
scalar field potential. 

We stress that the gravitating fermionic lumps described in the present paper 
exist when there appears a false vacuum and there are fermions coupled to the 
relevant scalar field. Recently it has been shown that such a false vacuum appears 
in the early universe in the standard Einstein- Weinberg-Salam theory of electroweak 
and gravitational interactions, if the universe has a spatial section as in the closed 
Friedmann-Robertson- Walker universe. In the early stage, nontrivial gauge fields 
yield a false vacuum in the Higgs field, to which quarks and leptons have Yukawa 
couplings. In other words gravitating fermionic lumps may be copiously produced 
in the framework of the standard model. 

In §2 we set up the problem and introduce useful variables in terms of which the 
field equations are written. In §3 we give the energetics of fermionic lumps, based 
mainly on flat spacetime. This illustrates how such lumps become possible when 
the scalar field potential has a false vacuum configuration. In §4 the behavior of the 
solutions is investigated analytically inside the lump, numerically in the transition 
region, and analytically outside the lump. More details of the solutions are given 
in §5 with a focus on the dependence of the solutions on various parameters of the 
model. It is seen how the solutions change as the energy scale of the model is lowered. 
A summary and conclusions are given in §6. 



We consider a real scalar field in the Einstein gravity whose Lagrangian is given 



where R and V[(j)] are the scalar curvature and the scalar potential, respectively. The 
last term in ()2-l|) represents a source for the scalar field. It naturally arises if there is 
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Fig. 1. A scalar field potential V[(j>\ with two minima. The quantity (j)ps = g(j>{tpil') at the core is 

also displayed for g^'ij)^) > 0. In the lump solutions described in the present paper, gl.fiKV'V') 
is much larger than the energy density of the false vacuum at /i . 



a fermion field ^ with a Yukawa interaction —gcpijjil)}^'^ In this case ps{x) = giipil^). 
As shown schematically in Fig. ^ ^[0] has two minima at /i and f2- We take 

V[4>] = ^(0 - /2)|<A' - ^(/2 + 4/l)0' - ^/2(/2 - 2/l)(0 + /2) 

V'[cP] = X(t>{^-h){^-f2). (2-2) 

In the potential V[(p], the false vacuum appears at (p = fi and the true vacuum at 
(j) = f2- We investigate spherically symmetric configurations in which the scalar field 
varies from the false vacuum to the true vacuum as the radius r increases. Solutions 
0(r) are to satisfy the two conditions (i) (/>(0) ~ /i and (ii) 4>{oo) = f2- 

In H2-l() we have added the interaction term —(j)pg{x) to the scalar field. This 
interaction is vital to have solutions satisfying the condition (ii). As we shall see 
below, solutions regular at the origin are uniquely specified by the value of 0(0) with 
a given ps{x). Without the coupling to a source —(j)ps{x), the configuration (j){r) 
either settles at = after oscillation as r — > oo, which corresponds to a totally 
unstable configuration, or comes back to (/> = /i, or overshoots 4> = fi to diverge to 
— oo as r approaches oo. Solutions of the second type indeed exist. They are cosmic 
shells.''^ The source term —(j)ps{x) gives a x-dependent linear potential for the scalar 
field. (See Fig. [H) 

We treat the system classically. If the source comprises fermions, ps{x) = 
g{i{ji{j{x) ) where g is a Yukawa coupling constant. If the fermions at the core are non- 

relativistic, {7pip{x) ) can be approximated by the fermion density {ip^ip{x) ) = Pn{x). 
In the present paper we suppose that the fermion density is small enough so that the 
back-reaction to the metric can be ignored. This condition is met in a simple model 
if the fermion couples to a pair of scalar fields in the specific manner described in the 
next section. A full treatment of the dynamics of the sources including back-reaction 
to the Einstein equations is left for a future study. 
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We look for spherically symmetric, static configurations. The metric of spacetime 
is written as 

ds^ = ^dt" - ^ - r\de^ + W 9dip^) (2-3) 
H 

where H and p are functions depending only on r. In the tetrad basis 



y^H 1 

eo = - — dt , ei = —=dr , 62 = rd6 , 63 = rs'mOdip , (2'4) 
P VH 

the energy-momentum tensors Tab = ea^e^^T^i, are 
Too = ^H^'^ + y[</>] + cPps , 

Tn = - V[(t>] - (fps , 

T22 = 133 = - n4>] - <^PS , 

others = , (2-5) 



where the prime represents a r-derivative. 
The Einstein equations are 



- — + ^^ = 87rGroo , (2-6) 

l(m^„.)_^^^S.aT,,, (2.7) 

r \ p J 

2 or \ p / r 



It follows from H2-6() that 



H = 1 , M = I Anr'^dr Tqq ■ (2-c 



r 

Adding (ES} to (EUl, we find 

„/ 

£1 = -4^Gr</.'^ . (2-10) 
P 

p{r) is a monotonically decreasing function. The equation of motion for the scalar 
field is 

Eq. (E31) follows from Eqs. (ESI, (Eill, and (EHH). 

After eliminating p{r), the equations are reduced to (|2-9p and 

</>"(r) + + 47rGr</.'(r)2 + 0'(r) = 1 (y' [0] + ^5) • (2-12) 
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In interpreting Eq. (|2-12|) it is convenient to introduce the new coordinate ^ as 

In the examples discussed below H{r) remains positive definite so that ^(r) is well 
defined. Then the equation reads 

^ + ne)^ = ^'M + P5(e) , (2-14) 

r(O = ^0 + 47rGr</.'(rf + . (2-15) 

In terms of the "time" coordinate ^ a particle with position moves in the po- 
tential — with a time-dependent external force ps{C)- There is friction F which 
can be either positive or negative. The solution we are looking for corresponds to 
a particle which starts from ~ /i and moves to /2 asymptotically. It must gain an 
energy as — < — V^[/2], which is possible as there is an additional force given 
by ps and F can be negative. 

Equations (|2-9() and (|2-12() define a set of nonlinear equations. To be definite we 
suppose that the source to (j) is given by 

ps{r) = PQe{Ri-r) . (2-16) 

We split the space into three regions: 

I. < r < i?i 
II. Ri<r<R2 
III. R2<r <(yo 

In region I, the spacetime is approximately anti-de Sitter with ps > and (po very 
close to, but still greater than, the location fs of the minimum of V[(j)] + (pps- 
(/,(0) = fs + S(j){0) with < 6(j){0)/\fs\ <C 1. It turns out that (j){r) varies little 
from fs in region I so that the equation of motion for cp may be linearized in this 
region. In region II, ps = 0. In this region the field varies significantly so that the 
full set of nonlinear equations must be solved numerically. In region III, ps = and 
the spacetime is approximately Schwarzschild. In this paper we focus on the case in 
which l/il ~ /2 ~ / = + /2), fr = (/2 - l/il)// « 1 and l/i - fs\ « / so that 

the linearization of Eq. (|2-12() is valid. The behavior of a solution (j){r) is displayed 
schematically in Fig. [3 



§3. Energetics of fermionic lumps 



It is instructive to see how fermionic lumps become possible from consideration 
of energetics in flat spacetime. The basic idea is that there reside nonrelativistic 
fermions at the core, which generate an additional linear interaction (pps- For non- 
relativistic fermions (ipip) ~ (ip^ip) = Pn so that the Yukawa interaction gcpipip 
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Fig. 2. The behavior of cj>{r). In the lump solutions \fi — fs\ ^ and R2 — Ri ^ -Ri. 
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Fig. 3. The behavior of H{r). H'{r) is discontinuous at _Ri as Too is discontinuous there when ps 
is given by 12-161 1. 



generates a linear potential gpn(l> for the scalar field inside the lump. Let R be the 
radius of the lump. We suppose that the total fermion number Np is conserved. Np 
is kept fixed in the consideration below. 

The scalar potential V^[</>] takes the form depicted in Fig. ^ The total potential 
V[(p] + gpn4' lia-s a minimum at fs < 0. It is supposed that V[fi] = eo > but the 
energy density at the minimum V[fs] + gpnfs = e < becomes negative for the 
lump solution, fs and e depends on pn or equivalently on R. 

Inside the lump r < R, the scalar field satisfies (j) ^ fs whereas outside the lump 
4> = f2- A degenerate nonrelativistic fermion gas has energy density £f = mc^pn + 
Apn^^ where A = 3{3TT^)^/^h'^ /lOm. The total energy of the lump is approximately 
given by 

E{R) = {Sf + e{pn)] ^ + 47ri?V + Ef^^.,n{R) (3-1) 

where a is the surface tension resulting from varying (p in the boundary wall region. 
-S/,wall(-R) is the contribution to the energy from fermions localized in the boundary 
wall region. It has been estimated in ref.^^^ to be about 2N'^J^^y^/2>R where Nf^^^^w 
is the number of fermions confined in the boundary wall region. 
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Fig. 4. The behavior of p(r). 



Let the total fermion number Np ~ (47ri?^/3)/9n + -^/,waii be fixed. e(p„) ~ eo 
for p„ ~ or as i? — oo. For large i?, ~ (47reo/3)i?^. For large p„ (small R) 

the total potential is approximated by (A/4)(/)^ + (7Pn0 near the minimum so that 
e ~ —\\~^^^{gpnY^^- When Np ~ (47ri?^/3)p„, the <Sj-term dominates over the 
e-term and E{R) ~ A'R''^ {A' > 0) for smah R. In the case that Pr, remains small 
for small i?, Ef^^gii{R) becomes important and E{R) ~ A"R~^ {A" > 0). In either 
case E{R) has a minimum, say, at R where R is the size of the lump. 

Of course the radius R cannot be too small in order for the nonrclativistic 
approximation to be valid. Further the above fermion lump configuration need to 
be energetically favored over a configuration in which fermions reside in the scalar 
configuration </> ~ /2- 

The above-stated conditions can be satisfied if the mass m of the fermion in Sf 
differs from the mass mo in the vacuum = /2. We suppose that e < inside the 
lump. In this situation it is energetically favorable for fermions to reside inside the 
lump rather than in the boundary region so that Np ~ (47ri?^/3)p„. For a Fermi 
ball discussed in refs.^^^'^^^ A''^;' = A^/,waii and p„ = 0. Note that e < implies that 

Pn > ^ - - ffr (3-2) 

91 9 

where / = \{\fi\ +/2) and = (/2 — |/i|)//- The nonrelativistic approximation for 
the fermi gas is justified if mc pn dominates over in £f; 

m»pV3>(M:)^/^/. (3-3) 

We would like to have the anti-de Sitter space inside the lump; + e < 0, or more 
strongly £f <C |e|. This leads to 

M~5/»m . (3-4) 

Pn 

In the examples of the lump solutions given below we shall have |e| eo- In other 
words gpn\fs\ S> Nevertheless fs is close to /i, and it is useful to introduce 
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the parameter hs = (/i — fs)/ f ^ 1- Pn and hs are related by gpn ~ 2A/^/i5. The 
conditions and (jTlj) read 

2hs>fr , i-^j • (3-5) 

It also implies that Pnl f'^ ~ '^^hg/g <^ 2g^. 

To satisfy (|3-5|) the fermion mass inside the lump, m, must be much smaller 

than gf. If the fermion acquired a mass only from the coupling gcpil^ij)^ then we 
would have m ~ gf . This apparent dilemma can be circumvented if the fermion 
has Yukawa couplings to more than one scalar field. Suppose there are two scalar 

fields, (j) and <P, which couple to the fermion through [gcj) + G^)'ilj'ip. In the true 

vacuum {(p) = f2 and (^) = /s. Hence the fermion mass in the true vacuum is 
rriQ = 5/2 + G/3 whereas inside the lump m = gfi + G/3 . We imagine that the two 
terms nearly cancel each other inside the lump such that m <^ gf and mo ~ 2gf. 
This scenario, at the same time, provides the stability of fermions inside the lump. 
It costs a huge amount of energy for a fermion to escape outside the lump. 

We stress that although we suppose, to facilitate numerical evaluation, that 
m ^ rriQ in the subsequent discussions, this condition is not absolutely necessary for 
the lump solutions to exist. In general cases we need to treat fermion contributions 
more accurately, taking into account the back-reaction to the metric as well. 



§4. Behavior of solutions 



Having made the assumption described in the previous section, we come back 
to the problem of solving the equations (|2-9|1 and 

(i) Region I 

Near the origin 0(r) is very close to fs which is given by 

V'[fs]+Po = . (4-1) 

fs itself is close to /i. Denoting 0(0) by (po, we have < (po — fs \fs\- The 
regularity of the solution at the origin leads to 

^{r) = <Po + <P2r^ + ■■■ , </>2 = ^ {V'[(l>o] + po) , 

p{r) = 1 + P4r^ H , p4 = -AttG(PI 

M{r) = m3r^ + ■■■ , ^3 = ^^(f[(/.o] + /Jq-^o) 

H{r) = 1 - IGmsr^ + ■■■ . (4-2) 
In region I the spacetime is approximately anti-de Sitter. As (/> ~ fs, 



roo = e = V[fs]+ pofs <0 , ^ = 1 + ^ ' "/ = V8^ ' P = ■ ^^'^^ 
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We have supposed that £f <^ \e\. The equation for (j){r) can be linearized in 5(f){r) 



(j){r) — fs- In terms of z = r'^ /a'j 



zil + z)^+(l + lz]^- \u;'a'f = , (4-4) 



d^^\2^rjd-z~r''f^ 

where = V"[fs]- This is Gauss' hyper geometric equation. The solution regular 
at r = is 

6cf>{r) = Sm . F(| + K, f - K, I; -z) , (4-5) 

where 



K=l^u;2a2+|. (4-6) 

We shall soon see that a solution with lump structure appears for uJUf ^ 1 with a 
particular choice of 6(j){0). The ratio of 5(j)'{r) to 5(f){r) is given by 

5(t>'{r) _ 4r / , 9 x F(| + | - |; -z) _ 2r 

6Hr) Sa^r 16)Fil + ^,l-n,l;-z)~ a}'^'^ ■ ^^''^ 

The deviation from fs at the origin, (5(/>(0), must be very small in order to have 
an acceptable solution. The behavior of the hyper geometric function for k 3> 1 and 
2; > is given 

r{ ^1 1 1 

F(a + K,a-K,c;-z) ~ — ^ K^^'^ z"i+3(l + 2;)t"3~'' exp{2Ksinh"^^| . (4-8) 

2v7r 

The ratio (5(/)(r)/(5(/)(0) grows exponentially as (4k)~^2;~^/^(1+z)~^/^ exp|2Ksinh~^ \/^}- 
Near r = R, 6(j)/\fs\ must be very small for the linearization to be valid. The ratio 
of F'{z) to F{z), J{z) in (jHH), is given by 

J{z) ^ for K > 1 , z>0 . (4-9) 

y/z{l + z) 



(ii) Region II 

In region II, </>(r) varies substantially and the nonlinearity of the equations plays 
an essential role. In this region the equations must be solved numerically. With fine 
tuning of the value of 5(f){Ri) nontrivial lump solutions are found. 

The algorithm for numerically finding solutions is the following. First (50(i?i) 
and fs are chosen and d(j)'{Ri) is evaluated using 1)4 -Tf) and H4-9p . In this paper we 
investigate the solutions in the case that — /s| <C / and the linearization in Region 
I is valid. To a good approximation the metric is given by H{Ri) = 1 + (i?i/a/)^ 
and p{Ri) = 1. With the boundary conditions 5(j){Ri) and 5(j)'{Ri) the equations 
are numerically solved in region II. The width of the transition region, w = R2 — Ri 
is approximately given by 1/\/A/-'^^ 
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S.OOxlo' S.Olxlo' 8.O2xl0' 

r 

Fig. 5. (t>{r) of a solution with //Mpi = 0.0002, /,. = 0.0002, A = 1, /is = - fs)/f = 0.005 and 
Ri/lpi = 8 X 10^. (/> and r are in units of J\fpi and Ipi, respectively. 



1.86r— r 




1.83- 



' 8.00xl0' 8.01xl0' 8.02xl0' 

r 

Fig. 6. H{r) of a solution with //A/pi = 0.0002, fr = 0.0002, A = 1, /is = - /s)// = 0.005 and 
Ri/lpi = 8 X lO'^. r is in units of Ipi. 

The behavior of solutions in region II is displayed in Fig. |3 When the values of 
the input parameters are chosen to be A = 1, f /M-p\ = 0.0002, fr = Af j f = 0.0002 
and hs = {fi — fs)/f = 0.005, then the output parameters are polpi = 8.0593 x 10^^^, 
e/M^j = -1.5944 x lO^^^ |e|/eo = 74.738, af/lpi = 8.6525 x lO'^ and k = 1.23275 x 
lOl Here hi is the Planck length. We note that po/f = 0.0101, e//^ = -0.00997, 
and ajf = 17305. For Ri/lpi = 8 x 10^ {Ri/af = 0.924589) we find a solution with 
6(p{Ri) / Mp\ = 5.545749- • • x 10~^. In this example, the value of 6(/) at the origin 
(r = 0) is found from Eq. ()4-5() to be 9.5 x 10^^^^^, which explains why numerical 
integration of (j){r) from r = to i?i is not possible. A small discontinuity in (f)" 
appears at r = Ri due to the discontinuous change in ps{r). 

The field (/>(r) approaches /2 for r > R2. In the numerical integration Ri and 
fs are kept fixed while 5(j){Ri) is varied. If (5</>(i?i) is taken to be slightly smaller, 
then 4){r) starts to deviate from /2 in the negative direction, approaching fi as r 
increases. If 5(l){Ri) is taken to be slightly larger, then (p[r) overshoots /2, diverging 
to +00 as r increases. With just the right value of 6(j){Ri), the spacetime becomes 
nearly flat in region III. 
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8.02x10 



Fig. 7. p(r) of a solution with f/Mm = 0.0002, fr = 0.0002, X = 1, hs = (/i - fs)/f = 0.005 and 
Ri/lpi = 8 X 10^. r is in units of Ipi. 



o 5x10" 




8.00x10 



8.02x10 



Fig. 8. The energy density Too for a solution with //Mpi = 0.0002, = 0.0002, A = 1, /ig = 



(/i-/s)// = 0.005 and i?i//pi 
For r < Ri, Too ~ e < 0. 



8x 10'. Too and r are in the units of Afpj and Ipi, respectively. 



(iii) Region III 

The behavior of the solution in region III is easily inferred. From the numerical 
integration in region II both H2 = H{R2) and p2 = p{R2) are determined. In region 
III the metric can be written in the Schwarzschild form 



H{r) = 1 
p{r) = p2 



2GM 



(4-10) 



where M is the mass of the lump. As (j){r) approaches its asymptotic value /2 and 
Too vanishes, H{r) must take this form. The value of M is determined numerically 
either by integrating Too in Eq. (|2-9() . or by fitting H{r) in Eq. H2-9() just outside the 
shell. In the numerical example presented above, M < 0. As can be seen in fig. El 
H{r) converges to the asymptotic value from above. A rough estimate of M is 
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8.00x10 8.01x10 8.02x10 

r 



Fig. 9. The radial pressure Tn for a solution with f/Mn = 0.0002, /r = 0.0002, A = 1, hs = 
(/i-/s)// = 0.005 andiii/Zpi = 8x10^ Tn and r are in the units of A/pj and Ipi, respectively. 
For r < Ri, Tn ~ -e > 0. 



§5. Scales of the solutions 

In the previous section we have solved the linearized field equation for 0(r) in 
region I where ps ^ and the deviation from fs was found to be small. We nu- 
merically determined the behavior in the nonlinear region II. The resulting structure 
is a lump with negative mass and fermions at its core. In this section we present 
more detailed numerical results. As discussed in the previous section the boundary 
between regions I and II, located at the matching radius Ri, is rather arbitrary, 
subject only to the condition that the linearization be accurate up to that radius. 
Precise tuning is necessary for Ri and (5(/)(i?i) in order to obtain a lump solution. 
Technically it is easier to keep Ri fixed and adjust 54> at Ri. If 64>{Ri) is chosen 
too small, (p comes back toward /i, but cannot reach it, eventually oscillating about 
= as r increases. If 5(p{Ri) is chosen too large, (p overshoots /2 and continues to 
increase. With just the right value of 54){Ri), (j) will approach /2. There can appear 
shell solutions in which (p goes back to /i at r = oo. Such solutions are discussed in 
ref.7) 

One example of the lump solutions is displayed in Fig. [S] for the parameter 
values //Mpi = 0.0002, fr = 0.0002, A = 1 and /i^ = {fi-fs)/f = 0.005. If we take 
Ri/lpi = 8 X 10^ {Ri/af = 0.924589), 6(j){Ri) musto be fine-tuned to more than 7 
digits: 6(p{Ri)/Mpi ~ 5.545749- •• x 10~^. In the transition region both H{r) and 
p{r) decrease in a one-step fashion. See Figs. iniand[7| Inside the lump H{r) is given 
by H4-3|l . whereas outside the lump it is given by 

2GM , , 

H = \ (5-1) 

r 

where M/Mpi k. -3.337 x lO'^ (M// k. -1.669 x 10^^). p(r) assumes the constant 
values 1 and 0.9944 inside and outside the lump, respectively. 

The behavior of the energy-momentum tensors, Tqo = —T22 = — 233 and Tn, are 
displayed in Fig. Eland Fig. 1^1 The energy density Tqo has one sharp peak associated 
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0.04 



0.025 h 

0.00016 0.00024 0.00032 0.0004 
/ 

Fig. 10. The / dependence of 5(i>{Rx). fr = 0.0002, Ri/af = 0.9 and hs = (/i - fs)/f = 0.005 are 
fixed. / is in a unit of Afpi. 



with the rapid variation of (j)- The radial pressure Tn is ~ — e > for r < Ri and 
becomes negative at Ri. It increases to zero quickly, and it remains zero outside 
the lump. The absolute value of Tii(i?i) in the transition region II is very small 
(~ 2.3 X IQ-^^M^j) compared with the maximum value of Tqq (~ 8 x 10 ^^Mpj). 
The contributions of the kinetic energy, ^Hcj)''^, and potential energy, V[(j)], almost 
cancel each other in the transition region. 

The model contains several dimensionless parameters, A, f/Mpi, fr, Ri/af and 
Pof^- When the values of these parameters are varied, the size of the resulting lump 
solutions also varies. In the numerical evaluation we took values of f /M-p\ and fr in 
the range between 10""^ and 10~^. It is of great interest to determine the structure 
of the lump solutions when, for instance, f /M-p\ = 10^^^. One can obtain insight 
into this problem by investigating the dependence of the solutions on the above 
parameters. 

In the numerical investigation values of A, f /Mp\, fr, Ri/aj and hs are given, 
and we try to find a desired value of 6(j){Ri) for a solution to exist. In Fig. IIUI 
6(j){Ri)/\fs\ is plotted as a function of /. It is seen that as / increases, the value 
of 5(p at Ri needs to be increased in order to get a lump solution. The numerical 
evaluation becomes unreliable when 5(j){Ri)/\fs\ becomes large and the linearization 
in Eq. I|2-12() is no longer valid. 

In Fig. [m 5(l){Ri)/\fs\ is plotted as a function of fr- We find that 6(j){Ri)/\fs\ 
increases as fr increases. In Fig. 6<j){Ri)/\fs\ is plotted versus Ri/aj. As Ri/aj 
increases, 54>{Ri)/\fs\ decreases and approaches a constant (~ 0.022). 

In Fig. d(j){Ri)/\fs\ is plotted versus hs- In this figure 6(j){Ri)/\fs\ decreases 
as hs increases and then it begins to increase when hs becomes greater than approx- 
imately 0.015. The reason for the increase in the value of 6(j){Ri)/\fs\ is that at some 
point 1/5 — fi\ becomes greater than 5(j){Ri), which means that (j){Ri) < fi- In this 
situation the solution will diverge to negative infinity if we do not take 6^{Ri) large 
enough. In fig. El -^i is plotted as a function of / with other parameters, including 
5,/.(i?i)/|/5|, kept fixed. 

With the above results we can obtain good information concerning the mass 
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0.0002 0.0004 0.0006 0.0008 



Fig. 11. The fr dependence of 5(j}{Ri). f = 0.0002, Ri/af = 0.9 and hs = (/i - fs)/f = 0.005 are 
fixed, fr is in a unit of Mpi. 




4 



Fig. 12. The _Ri dependence ofS(j){Ri). f/Mpi = 0.0002, fr = 0.0002 and hs = (/i -/s)// = 0.005 
are fixed. 

M of the lumps. In Figs. El and M is plotted as a function of / and 
respectively. We find that the estimate H4-lip is good. As an example, let us take 
fr = 0.0002, hs = 0.005, A = 0.01, and g = I. Then a/ ~ 1.2 • 10^ • (lGeV//)2m 
and -M ~ 4.6 • 10^^ • (IGeV//)^ (i?i/o/)3kg. As / becomes small, the size of the 
lump becomes very large. 

§6. Discussion 

In this paper we have demonstrated that there exist fermionic gravitating lumps 
with a false vacuum core of the scalar field. This is a curious structure which may 
appear at various scales. We have constructed lump solutions for the case hs = 
(/i - fs)/f >fr = (/2 - l/il)//, which implies that Xfr/g < Pn/2f « g\ (See Eq. 
H3-5() .l If the typical energy scale is very high, the size of the lump will be small but 
it can be produced abundantly in the early universe. If the energy scale becomes 
low, this size may increases to a cosmic scale. 

We have focused on the cases in which fermion contributions or back-reaction 
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0.005 0.01 0.015 0.02 0.025 

Fig. 13. The hs dependence of 5(t){Ri). f/Mm = 0.0002, fr = 0.0002 and Ri/aj = 0.9 are fixed. 

1.1 
1.0 
0.9 

°= 0.8 
0.7 

0.00018 0.0002 0.00022 0.00024 0.00026 
/ 

Fig. 14. Ri versus /. fr = 0.0002, hs = (/i - fs)/f = 0.005 and 5<^(i?i)/|/s| = 0.03 are fixed. / 
is in a unit of Afpi. 

to the metric is small. Certainly this restriction needs to be relaxed in order to 
accommodate more realistic models in the particle theory. The lump solutions pre- 
sented in this paper have negative masses as seen outside the lumps. As the mass 
and density of fermions inside the lump increase, the mass of the lump will become 
positive. The solution probably continues to exist in this case, though explicit con- 
struction of such a solution is necessary by taking account of the back-reaction to 
the energy-momentum tensors. An investigation along this line is in progress, and 
we hope to report it in the near future. 

A false vacuum state may appear in the standard model of electroweak and 
gravitational interactions. In ref. 17) it is shown that the false vacuum of the Higgs 
field emerges in the early universe as a result of the winding of the gauge fields in the 
X Robertson- Walker spacetime. As quarks and leptons have relevant Yukawa 
couplings to the Higgs field, gravitating lumps with quarks and leptons at the core 
may be copiously produced. As the universe expands, the barrier separating the false 
vacuum from the true vacuum disappears so that the gravitating fermionic lumps 
would become unstable, the fermions inside the lumps dissipating to infinity. It is of 
great interest to elucidate the consequences of this process. 
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Fig. 15. M versus /. fr = 0.0002, hs = (,fi - fs)/.f = 0.005 and Ri/af = 0.9 are fixed. / and M 
are in a unit of Mpi. 
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0.6 0.7 0.8 0.9 1 

Fig. 16. M versus Ri. //Mpi = 0.0002, fr = 0.0002 and hs = 0.005 are fixed. M is in a unit of 
Mpi. 

Furthermore, in the higher-dimensional gauge theory defined on orbifolds, false 
vacua naturally appear in the gauge field configurations. ^'^^ It would be interesting 
to investigate if the gauge interactions of fermions produce gravitating lumps, as we 
have found for the Yukawa interactions. 
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